Frustration-induced insulating chiral spin state in itinerant triangular-lattice magnets 
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We study the double-exchange model at half-filling with competing superexchange interactions 
on a triangular lattice, combining exact diagonalization and Monte-Carlo methods. We find that in 
between the expected itinerant ferromagnetic and 120° Yafet-Kittel phases a robust scalar-chiral, 
insulating spin state emerges. At finite temperatures the ferromagnet - scalar-chiral quantum critical 
point is characterized by anomalous bad-metal behavior in charge transport as observed in frustrated 
itinerant magnets R2M02O7. 



Geometric frustration is encountered in numerous 
magnetic condensed matter systems and gives rise to a 
variety of fascinating magnetic phases, such as spin-glass, 
spin- liquid, and spin-ice [l|, Q- Traditionally geometri- 
cally frustrated magnetism is studied in the context of 
wide bandgap materials with localized spins. Recently 
the interest broadened to frustrated magnets with metal- 
lic character such as Tl2Mn207 and R2M02O7 (R: rare 
earth ion) [j-Q- In these metallic systems, frustration in 
the magnetic sector has a strong impact on charge dy- 
namics, often resulting in heavy-fermion type behavior 
in transport J7l-[l0|. 

Besides being of geometric origin, magnetic frustra- 
tions can also be due to a direct competition between 
ferromagnetic (FM) and antiferromagnetic (AFM) inter- 
actions. Perhaps the most elementary example is FM 
double- exchange (DE) competing with AFM superex- 
change (SE) [11[. Magnetic competition of this kind has 
been studied extensively on hyper-cubic lattices in the 
context of the colossal magnetoresistance manganites, 
where depending on carrier concentration the competing 
interactions can induce new forms of magnetic order or 
phase-separation |12h14| |. Competing FM and AFM in- 
teractions can in principle also lead to non-collinear and 
even non-coplanar magnetic states, which are raising in- 
terest in the rapidly emerging and seemingly disparate 
fields of multiferroics and topological insulators |15l - [l7| . 

In an interesting class of materials both types of frus- 
trations - geometric frustration and competing FM- and 
AFM interactions - are present and give rise to a set of in- 
triguing physical properties [ISl |l9[ . Recent experiments 
on pyrochlore R2M02O7, for instance, show a transition 
from a FM metallic to a spin glass insulating state as R 
changes from Nd to Dy Q. This transition can be con- 
trolled by external pressure. In the vicinity of the transi- 
tion an unusual diffusive metallic state appears, showing 
a temperature independent resistivity down to very low 
temperatures [8| . A similar competition between FM and 
AFM interactions is also relevant in the triangular lattice 
systems such as Gdl2 and H-doped Gdl2 [2(]| - |23j . 

In this letter, we study in a wide temperature range 
the double-exchange (DE) model on a triangular lat- 
tice in the presence of frustrating AFM SE interactions. 
Three magnetically ordered groundstates are present in 



the phase diagram: an itinerant FM state on one end, the 
120° Yafet-Kittel phase on the other and an insulating 
(I) non-coplanar scalar chiral (SC) state in between 24 1. 
Apparently the SC-I state emerges from the competition 
between the two metalhc phases: FM and 120°. The FM 
and SC state are separated by a quantum critical point 
(QCP), whereas the transition from the SC to the 120° 
phase is first order in nature. At finite temperature in 
the vicinity of the QCP the system is characterized by a 
temperature-independent resistivity and a linear inverse 
magnetic susceptibility, down to very low temperatures. 
This bad metallic behavior closely resembles the experi- 
mental observations on the R2M02O7 compounds. 

We consider the elementary one-band DE Hamiltonian 
in the presence of antiferromagnetic SE interactions, in- 
troduced by De Gennes in the 1960's [Il[ - but now on 
a frustrated triangular lattice. The full Hamiltonian is 



H = -J2*^J{4c,+H.C. 
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where Cj and c] are annihilation and creation operators 
for electrons with spin parallel to the core spin Si. (ij) 
denotes the nearest neighbor (nn) pairs of sites on a 
triangular lattice. Js denotes the strength of AF cou- 
pling between nn core spins, tij = t[cos{6i/2) cos{9j/2) + 
sin(6'j/2)sin(6'j/2)e-'('^'-'^^^] denote the hopping ampli- 
tudes which depend on the polar and azimuthal angles 
{9i,4'i,9j,(j)j} of the nn core spins due to the double- 
exchange mechanism. All energies are measured in units 
of the hopping parameter t. The core spins are classi- 
cal unit vectors and we focus on the case of half filling. 
The model Eq. (1) can be obtained as the large cou- 
pling limit of the Kondo lattice model (KLM) , Hxlm = 

-iE(,,),.(cls- + ^-c-) + ■^K ^^ S* • ^^- The 1st term 
in Eq. (1) corresponds to \Jk\ — > 00 in Hklm, and the 
lowest order correction (^ i^/| J^l) leads to the superex- 
change term. Furthermore, the Hubbard model in the 
mean-field approximation has a similar structure as the 
KLM with the local moments originating self-consistently 
within a single band [l6[. Therefore, the results discussed 
in this paper also bear relevance for the related models. 
The only unbiased method to study this model is the 
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FIG. 1: (Color online) (a) Magnetization and (b) inverse mag- 
netic susceptibility as a function of temperature for various 
values of superexchange coupling Js- The inset in (b) shows 
the Curie- Weiss scale extracted from the x~ {T) as a func- 
tion of Js- (c) Spin-structure factor S'(q) at T = 0.001 along 
the principal symmetry directions in the momentum space for 
different values of Js- The inset shows the Brillouin zone for 
a triangular lattice and the path along which S(q) is plotted. 



hybrid scheme involving exact-diagonalization (ED) of 
the fermion problem and Monte Carlo (MC) for clas- 
sical core spins [12[. In order to achieve larger lattice 
sizes, which are essential for computing transport prop- 
erties, we employ the traveling cluster approximation 
(TCA) [25|. This method is based on the observation 
that changes in total energy induced by a change in clas- 
sical variable on a site can be estimated accurately by ED 
of a cluster Hamiltonian around that site. The method 
has been benchmarked for similar models, and has proved 
successful in the study of manganites [23, |2S| • Most of the 
results in this work are obtained on lattices with N = 24^ 
sites using a cluster size iVc = 6^. Typically ~ 10^ MC 
steps arc used for equilibration and a similar number of 
steps for computing thermal averages on classical spin 
variables. For electronic properties, which still requires 
the ED of the full Hamiltonian, ~ 10"^ steps are used. 

Fig. 1(a) shows the calculated temperature depen- 
dence of magnetization (m) = {\i/NJ2i^i)\)j where 
(...) denotes thermal averaging. In the absence of Js 
wc find a FM groundstate with a Curie temperature 
Tc = 0.15 as inferred from the inflection point in the 
{m)(T) curve. Upon increasing Js a monotonic reduc- 
tion in Tc is observed while the groundstate continues 
to be a saturated FM. For Js = 0.08 the ferromagnetism 
is destroyed by the competing AFM interactions and wc 
do not find any long-range order in the groundstate. In 
Fig. 1(b) we show the inverse of magnetic susceptibility 



computed as x = T~^{{m'^) — (m)^). The inverse suscep- 
tibility x~^ clearly follows a Curie- Weiss (CW) behavior 
(x oc 1/(T— 9)) above a characteristic temperature. The 
CW scale Q, which is obtained by extrapolating the lin- 
ear high-temperature behavior of x~^, is shown in the 
inset in Fig. 1(b). Q reduces monotonically and changes 
sign upon increasing Js, indicating a change in the na- 
ture of effective magnetic coupling from FM to AFM. A 
remarkable similarity of these results with those obtained 
on a pyrochlore lattice indicates that the essential physics 
is determined by the frustrated nature of the lattice [l^ . 

In order to identify non-trivial long-range ordered 
magnetic phases wc compute the spin structure factor, 
5(q) = l/N^J2{^i ■ Sj) e'lt'-'-'-^), where r,, Tj are the 
real-space location of spins Si, Sj . We show the ground- 
state structure factor along the symmetry directions in 
the momentum space in Fig. 1(c). For Jg = 0.04, the 
S'(q) has a single peak at q = (0, 0) as expected for a FM 
state. In the regime 6^0 {Js = 0.08), the structure 
factor has no peak at any q suggesting the absence of any 
long-range magnetic order. For Js = 0.16, we find peaks 
at q = (tt, 7r/-\/3) and (0, 27r/-\/3). This corresponds to a 
four-sublattice ordered state with a non-vanishing scalar 
chirality n = X](^i ' ^i ^ ^k), where the sum is over 
indices forming a triangle taken in the counterclockwise 
order. The onset temperature Tgc for the SC state can 
be inferred from the temperature-dependence of k. This 
state has been discussed for the KLM within mean-field 
and variational schemes, and has also been studied in 
the context of high-Tc su per conductors and recently the 
topological insulators [Ifl, [l3, [23 ■ Here we show that the 
SC state is the groundstate for the model Hamiltonian 
Eq. (1) over a window of phase space, and does not rely 
on Fermi surface nesting [131 ■ Finally as the Js becomes 
large the system should approach towards the classical 
120° state. This is seen in the S{q) data for Js = 0.44, 
where the peak is located at q = (47r/3, 0). 

To investigate the repercussions of the competing mag- 
netic interactions on electronic properties, we first fo- 
cus on the evolution with Js of the electronic density of 
states (DOS), D{E) = (Ek^i^ " ^k)), where we ap- 
proximate the (5-function by a Lorentzian of width 7: 
t:5{E - Efc) - 7/(72 + (£; _ gj^)2]^ ^he DOS for differ- 
ent temperatures and values of Js is shown in Fig. 2. 
For small Js the DOS does not show any drastic changes 
in shape upon reducing T (see Fig. 2(a)) and an increase 
in bandwidth is found as expected in a DE model. For 
Js = 0.08 the DOS at low temperatures is identical to 
that in the high temperature paramagnetic (FM) state 
(see Fig. 2(b)). This unusual behavior of the DOS is con- 
sistent with the facts that (i) x~^ follows the Curie- Weiss 
behavior down to T = 0, and (ii) the spin-state remains 
disordered down to low temperatures for Js = 0.08 (see 
Fig. 1). For Js = 0.16, where the groundstate has the 
peculiar non-coplanar order, an energy gap appears in 
the DOS near Ep at low temperatures (see Fig. 2(c)). 
This is consistent with the previous discussions on this 
non-trivial state [lajllJillJ- The effect of this energy gap 
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FIG. 2: (Color online) Electronic density of states for varying 
temperatures with (a) Js = 0.04, (b) Js = 0.08, (c) Js = 0.16 
and (d) Js ~ 0.36. A Lorentzian broadening 7 = 0.04 is used. 



persists at higher temperatures in the form of a pseudo- 
gap feature at Ep. For larger Js an unusual flip in the 
DOS shape w.r.t. the high temperature DOS occurs near 
T = 0.03 (sec Fig. 2(d)) [H]. We find that this flip in 
the shape of DOS is correlated with the rise in S'(q) at 
q = (47r/3, 0) and hence indicates a 120° state. 

We compute the dc conductivity {a{u!)\uj^o) using the 
Kubo-Greenwood formula and the exact eigenspectrum. 
In this work we use cr(wo)Uo=o.03 as an approximation 
to dc conductivity except for low-temperature metallic 
states, where we make use of the conductivity sum rule 
[29I |30 | . A typical double-exchange behavior in resistiv- 
ity is found in cases where the groundstate is FM (see 
Fig. 3(a)). For Js = 0.08, which corresponds to a dis- 
ordered magnetic groundstate as seen in S'(q), we flnd 
an unusually flat resistivity. Further increase in Js leads 
to an insulating behavior {dp/dT < 0) at low temper- 
atures with a diverging resistivity at T ^ 0.01. For 
Js = 0.30 one finds a sharp reduction in p at a char- 
acteristic temperature associated with the rise in S'(q) at 
q = (47r/3,0). dp/dT remains positive in the low tem- 
perature regime until the system enters the scalar chi- 
ral insulating state near T = 0.005. Eventually, beyond 
Js = 0.40 the groundstate becomes the 120°-state. 

In Fig. 3(b) we show resistivity as a function of Js for 
low and high temperatures. This data can be compared 
with the pressure dependence of resistivity reported in 
[8|, since application of external pressure alter the Js/t 
ratio and tunes the system across the FM to PM tran- 
sition. We find the resistivity at high temperature in- 
creases monotonically, whereas at low temperatures it 
shows a non-monotonic behavior. These features quali- 
tatively agree with the experimental observation [g . Wc 




FIG. 3: (Color online) (a) Resistivity as a function of temper- 
ature for various values of Js- (b) Js-dependence of resistivity 
p (filled symbols) and scaled DOS (open symbols) at Fermi 
level Ef for low and high temperatures. 



also plot the DOS at Ep as a function of Js for low and 
high temperatures. Naively at low temperatures one ex- 
pects the resistivity to be inversely related to the DOS 
at Ep, since only the states near Fermi level contribute. 
We indeed find this behavior at low temperatures: the 
peak is p corresponds to a dip in D{Ep). 

In Fig. 4(a) we plot the peak values in S'(q) at var- 
ious q at T = 0.005 as a function of Js- Peak in S'(q) 
at q = (0,0) is a measure for ferromagnetism, that at 
q = (47r/3, 0) corresponds to 120° state and simultane- 
ous peaks at q = (7r,7r/-\/3) and (0,27r/-\/3) indicates a 
SC state. Upon increasing Js, S'(0,0) vanishes before a 
rise in S{'K,^l/^/i) and 5(0, 27r/-\/3), indicating that the 
associated transition at Js '^ 0.07 is quantum critical in 
nature. On the other hand the transition between chiral 
and 120° state is 1st order. The dashed curve shows the 
Jg-dcpcndence of the scalar chirality k. The results are 
summarized in a phase diagram in Fig. 4(b). The ferro- 
magnetic transition temperature {Tc) is inferred from the 
(to) (r) , the onset of pseudogap behavior is determined 
by directly looking at the DOS at various temperatures, 
and the transition to the 120° state is determined from 
the T-dependcncc of the S'(q) at q = (47r/3, 0). The ab- 
solute value |6| of the Curie- Weiss temperature is plot- 
ted as an estimate of the effective magnetic interaction in 
the system. The FM scale closely follows the value of O, 
indicating the unfrustrated nature of the system in this 
regime. In the intermediate- Js regime the groundstate is 
SC-I and a pseudogap behavior appears just above Tgc- 
This simple looking model provides an elegant example 
where a competition between two metallic phases leads 
to an intermediate insulating phase. 

For a simple physical picture of the various magnetic 
states we derive an effective model for the DE term as 
-^e// = ^Y.{ij) KijCos[a.ij /2), with cos{aij) = S^ • Sj, 



and Kj^ 
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here the angular bracket de- 



notes the quantum expectation value in the groundstate 
[3l| . The phase $y in ti-j = |iij|e'*'J is trivial for the 
FM and 120° states, which leads to real-valued hopping 
parameters. Minimizing Jscosiaij) — KijCos(aij/2) for 




FIG. 4: (Color online) (a) Low-temperature spin structure 
factor ^(q) as a function of superexchange coupling Js for 
various q. The dashed line is the scalar chirality k. (b) The 
phase diagram of the model in the temperature- Js plane. The 
spin directions within a unit cell are schematically shown for 
the 120° state, and are indicated as un-normalized {Sx, Sy, Sz) 
on each site for the non-coplanar SC state. 



a single bond with Kij — KpM ^ 0.32t wc find that the 
FM state becomes unstable for Js > 0.25Kfm = 0.08i, 
which is very close to the computed transition point. In 



contrast, the SC-I state has non-trivial $y, which give 
rise to local magnetic field terms in the effective Hamilto- 
nian. The total phase around a closed plaquette forming 
the 4-site unit cell is tt, and a staggered sign of this flux 
between neighboring plaquettes explains the opening of 
energy gap in the spectrum of the SC state. The self- 
consistent emergence of this magnetic flux clearly relies 
on the presence of itinerant fcrmions. A similar chiral 
spin state has also been recently reported in a KLM on a 
pyrochore lattice [32] . We also find non-trivial spin states 
for the present model on a checkerboard lattice [33, in- 
dicating that the results discussed here in detail for the 
triangular lattice have broad implications. 

To conclude, a scalar-chiral insulating state emerges 
when double exchange and superexchange interactions 
compete on a triangular lattice. The transition from fer- 
romagnetic metal to a scalar-chiral insulator is quantum 
critical in nature with both Tc and Tgc vanishing at a 
critical value of Js- We find that the chiral order is 
a consequence of a self-consistently generated magnetic 
flux, which relies crucially on the presence of itinerant 
fcrmions. There is increasing experimental evidence that 
such chiral spin states lead to anomalous Hall effect via 
the hopping generated magnetic flelds [^, |3J| . The tran- 
sition from a FM-metal to a spin-glass insulating phase 
reported in the experiments on R2M02O7 is qualitatively 
similar to that from FM-M to the SC-I phase d]. The 
observed flat-resistivity in the vicinity of the transition 
is very well captured within our calculations. 

SK acknowledges stimulating discussions with A. 
Mukherjee, G. V. Pai, S. R. Hassan and A. Paramekanti. 
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